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ABSTRACT 

We present high-speed optical photometry of the anomalous X-ray pulsar IE 1048.1- 
5937 obtained with ULTRACAM on the 8.2-m Very Large Telescope in June 2007. 
We dete ct IE 1048.1-5937 a t a magnitude of i' — 25.3 ±0.2, consistent with the values 
found bv lWang et alj (|2008f ) and hence confirming their conclusion that the source was 
approximately 1 mag brighter than in 2003-2006 due to an on-going X-ray flare that 
started in March 2007. The increased source brightness enabled us to detect optical 
pulsations with an identical period (6.458 s) to the X-ray pulsations. The rms pulsed 
fraction in our data is 21 ± 7%, approximately the same as the 2-10 keV X-ray rms 
pulsed fraction. The optical and X-ray pulse profiles show similar morphologies and 
appear to be approximately in phase with each other, the latter lagging the former 
by only 0.06 ± 0.02 cycles. The optical pulsations in IE 1048.1-5937 are very similar 
in nature to those observed in 4U 0142+61. The implications of our observations for 
models of anomalous X-ray pulsars are discussed. 
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1 INTRODUCTION 

The anomalous X-ray pulsars (AXPs) are a small grourj]] 
of isolated neutron stars in which the X-ray luminosity far 
exceeds the energy available from the spin-down. The AXPs 
are generally believed to be magnetars, in which the ex- 
cess luminosity is powered by the dec ay of an ultra-strong 
magn etic field, in excess of 10 14 G ( see I Woods fc Thompson! 
120061 ). An alternative explanation is the fallback disc sce- 
nario, in which some of the supernova ejecta fails to escape 
and forms an accretion disc around the neutron star, pro- 
vidin g an extra source of energy to power the X-ray emis- 
sion jvan Paradiis et al.lll995l : IChatteriee et alj|2000t lAlparl 
120011) . 

One way of discriminating between the magnetar and 
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fallback disc models is via optical observations. The mag- 
netar model predicts any optical emission must be non- 
thermal and magnetospheric in origin. Four plausible mecha- 
nisms have been considered - coherent plasma emission, syn- 
chrotron emission from electrons with high Lorentz factors, 
cyclotron emission from ions in the outer magnetosphere 
and curvature emission from bun ched electron-positron pairs 
in the inner magnetosphere fsee iBeloborodov fc Thompson! 
(2007) and references therein). The fallback disc model, on 
the other hand, predicts any optical emission is produced by 
reprocessing of the X-r ay light in the dis c and/or thermal 
emission from the disc l|Perna et al.ll2000h . 

The first AXP to be de tected in the optical p art of the 



spectrum was 4U 0142+61 (|Hulleman et al 



2000 ). Optical 



for an up-to-date catalogue of all known AXPs, including the 
various wavelengths at which each has been detected. 



pulsat ions were discovered in 4U 0142+61 bv lKern fc Martini 
(2002), and the fact that these pulsations have the same pe- 
riod, morphology and phase as the X-rays, but with 5-7 
times greater pulsed fraction, was reported bv lDhillon et al.l 
(2005). These results provided strong support for the magne- 
tar model - pulsed optical emission is indicative of a mag- 
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Figure 1. Left: Summed i'-band image of the field around 
IE 1048.1-5937, with a total exposure time of 10684 s. For clar- 
ity, only a portion of one of the two ULTRACAM windows is 
shown. The positions of IE 1048.1— 5937 and the comparison star 
are indicated by circles near the centre and bottom of the image, 
respectively. The central box shows the portion of the field that 
is plotted at a higher contrast on the right. The orientation of the 
image is marked on the upper right-hand side. The pixel scale is 
0.156 arcseconds/pixel, hence the field of view in this image is 
36 X 30 arcseconds. The vertical banding is due to residual bias 
structure. Right: Higher contrast plot of a 7.5 X 7.5 arcsecond field 
around IE 1048.1— 5937, highlighting the detection of the pulsar 
in the i'-band 



netospheric origin and disc reprocessing is unlikely in this 
case, as the optical pulsed fraction is higher than the X-ray 
pulsed fraction and there is no time delay between the two. 
Although it is possible to contrive ways in which the fallback 
disc model is consistent with the optical pulsations observed 
in 4U 0142+61, e.g. by assuming that the X-ray pulse pro- 
file that we observe is different to the X-ray radiation seen 
by the disc due to orientation or beaming effec ts, or by in- 
voking a hybrid disc-magnetosphere model (see lErtan et al.l 
l|2007h and references therein), the weight of evidence from 
the optical and other wavel engths lies heavily on the side 
of the magnetar model (see iMereghettil (|2008h for a recent 
review) . 

The detection of optical pulsations in other AXPs would 
provide valuable confirmation, or otherwise, of the results 
for 4U 0142+61 discussed above. Only one other AXP has 
been unambiguously identified in the optical: IE 1048.1- 
5937 (jDurant fc van Kerkwiikl 120051 ) 1 . In this paper we re- 
port on the first detection of optical pulsations from AXP 
IE 1048.1-5937, obtained only ~3 months after a bright X- 
ray flare in March 2007. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of IE 1048.1-5937 p resented in this pa - 
per were obtained with ULTRACAM (|Dhillon et all 120071 ') 
at the Nasmyth focus of Melipal, the 8.2-m Unit 3 of the 
Very Large Telescope (VLT) in Chile. ULTRACAM is a 
CCD camera designed to provide imaging photometry at 
high temporal resolution in three different colours simulta- 
neously. The instrument provides a 2.66 arcminute field on 
its three 1024 x 1024 E2V 47-20 CCDs (i.e. 0.156 arcsec- 
onds/pixel). Incident light is first collimated and then split 
into three different beams using a pair of dichroic beamsplit- 



Table 1. X-ray ephemeris for IE 1048.1-5937 llDib et alj|2008h . 
The epoch of the frequency and frequency derivative measure- 
ments given below falls on the same night as our VLT observations 
(09/06/2007 = MJD 54 260). BMJD refers to the Baryccntric- 
corrected Modified Julian Date on the Barycentric Dynamical 
Timescale (TDB). The errors on the last two digits of each param- 
eter are given in parentheses. This ephemeris is valid for BMJD 
54 229.0 - 54 280.0. 



v (Hz) 0.1548479469(42) 

v (10^ 13 Hz s" 1 ) -5.413(53) 

Epoch (BMJD) 54 260.0 



ters. For the observations presented here, one beam was ded- 
icated to the SDSS u (\ e ff = 3543A) filter, another to the 
SDSS g (4770A) filter and the third to the SDSS (7625A) i' 
filter. Because ULTRACAM employs frame-transfer chips, 
the dead-time between exposures is negligible: we used UL- 
TRACAM in its two- windowed mode, each of 250 x 200 pix- 
els, resulting in an exposure time of 0.963 s and a dead- 
time of 0.024 s. A total of 11095 frames of IE 1048.1-5937 
were obtained on the night of 2007 June 9, with each frame 
time-stamped to a relative (i.e. frame-to-frame) accuracy of 
' 50 £ts and to an absolu te accuracy of ~ 1 m s using a ded- 



icated GPS system (see iDhillon et al.l l2007h. Ob servations 
of the SDSS standard G163-51 (|Smith et al-lfeool ) were also 



obtained to flux calibrate the data. The night was photomet- 
ric, with no moon and i'-band seeing of 0.65 arcseconds. The 
sum of the 11 095 frames in the i'-band is shown in figure [1] 
which can be compared to the finding charts presented by 
iDurant fe van Kerkwiikl (|2005l ). 

The data were reduced u sing the ULTRACAM pipeline 
software (|Dhillon et al.ll2007h . All frames were first debiased 
and then flat-fielded, the latter using the median of twilight 
sky frames taken with the telescope spiralling. Adopting the 
same success ful approach that w e used in our study of AXP 
4U 0142+61 |Dhillon et al.ll2005t ), we extracted light curves 
of IE 1048.1-5937 using two different techniques: 



2.1 Technique (i) 

As part of a long-term monitoring project, IE 1048.1- 
5937 has been observed regularly (up to three times per 
week) since 1997 with the Proportional Counter Array 
(PCA) on board the Rossi X-ray Timing Explorer (RXTE ) 
[Kaspi et all l200ll : iGavriil fc Kaspl |2004| ; iDib et al.|[200^ ]. 
The X-ray spin frequency and frequency derivative of 
IE 1048. 1-5937 for the night of our VLT observations (MJD 
54 260) are given in table 1. For the first light-curve ex- 
traction technique, we shifted and added each of the 11 095 
ULTRACAM frames into 10 evenly-spaced phase bins us- 
ing the epoch and spin frequency given in table 1, result- 
ing in 10 high signa l-to-noise dat a frames. An optimal pho- 
tometry algorithm (|Navlorlll998h was then used to extract 
the counts from IE 1048.1-5937 and an i' ~ 17 comparison 
star ~ 12 arcseconds to the south-east of the AXP (see 
figure [TJ , the latter acting as the reference for the profile 
fits and transparency-variation correction. The position of 
IE 1048.1-5937 relative to the comparison star was deter- 
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mined from a sum of all the images, and this offset was 
then held fixed during the reduction so as to avoid aperture 
centroiding problems. The sky level was determined from a 
clipped mean of the counts in an annulus surrounding the 
target stars and subtracted from the object counts. 

2.2 Technique (ii) 

The second approach we took to light curve extraction was 
identical to that described above, except we omitted the 
phase-binning step and simply performed optimal photom- 
etry on the 11 095 individual ULTRACAM data frames fol- 
lowed by a periodogram analysis of the resulting time series. 
In other words, we made no assumption about the spin pe- 
riod of IE 1048.1-5937. 

3 RESULTS 

3.1 Magnitudes 

We were unable to detect IE 1048.1-5937 in u and g', 
at a 3<t detection limit of it' > 25.7 and g' > 27.6, re- 
spectively. This is unsurprising gi ven the high visual ex- 
tinction to the object [Av = 4.9; iDurant fe van Kerkwiikl 
12006ft . We did, however, clearly detect IE 1048.1-5937 in 
i' at a mag nitude of i! = 25.3 ± 0.2, as shown in fig- 
urefTl lDurant fc van Kerkwiikl l|2005h found that IE 1048.1- 
5937 was at a magnitude of I = 26.2 ± 0.4 (i' ~ 26.10 on 
06/06/2003, just over one y ear after its first X -ray flare 
was observed in April 2002. IWang et alJ (|2008l ) observed 
IE 1048.1-5937 at a magnitude of I = 24.9 ±0.2 (i'~24.8) 2 
on 07/05/2007, just over one month after its second X-ray 
flare began on 2007 March 21. These authors also measured 
a limit of i' > 24.5 on 15/07/2007. Our i'-band magnitude, 
which is slightly fainter and was obtained slightly l ater than 
the first /-band measurement of IWang et al. l|2008h . implies 
that we observed IE 1048.1-5937 whilst it was still relatively 
bright, but decli ning, from the most recent X-ray flare (see 
iTam et al.ll2008h . 

3.2 Pulse profiles 

The two data reduction techniques described in section [2] 
result in two different pulse profiles for IE 1048.1-5937. 

3.2.1 Technique (i) 

The first technique produced the pulse profile shown by 
the solid line in the top panel of figure [2] The pulse pro- 
file exhibits a broad, single-humped structure with a peak 
around phase 0.5 and a minimum around phase 0. There 
is a great deal of similarity in the morphologies of the op- 
tical pulse profile shown in the top panel of figure [2] and 
the 2-10 keV X-ray pulse profile shown below it, where the 

2 Noting that the visual extinctions to IE 1048.1— 
5937 and 4U 0142+61 a re approximately the same 
iDurant fc van Kerkwiikl2006h . we have assumed that the colours 
of IE 1048.1-5937 are the same as 4U 0142+61 l lHulleman et al.l 
2004 iDhillon et all [2005) and then used the equations of 
Smith et al.l (120021 1 to convert from I to i' . 
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Figure 2. Top: The solid and dotted lines show optical pulse pro- 
files of IE 1048.1-5937 in the i'-band obtained using techniques 
(i) and (ii), respectively (see sections 12.11 and 12,20 , Each pulse 
profile was first corrected for transparency variations using the 
comparison star shown in figure [T] although the correction made 
only a negligible difference to the light curves. The pulse profiles 
were then normalised by dividing by the mean number of counts. 
For clarity, two cycles are shown. Bottom: Averaged X-ray pulse 
profile of IE 1048.1— 5937 in the 2-10 keV energy ban d spanning 
the epoch of the optical observations l|Dib et al]|2008l) . Note that 
it is not possible to estimate the X-ray pulsed fraction from this 
profile as the PCA on RXTE has a 1° field of view and no imag- 
ing capability, rendering the background level uncertain. For this 
reason, no scale is given on the ordinate. 



latter is the average of the X-ray light curves in the pe- 
riod 23/05/2007-21/06/2007 obtai ned as par t of th e RXTE 
monitoring campaign described bv lDib et al] l)2008h . with a 
total effective integration time of 59.5 ks. Both profiles show 
the same broad, single-humped morphology. Moreover, since 
the X-ray light curve shown in figure[2]has also been phased 
using the ephemeris given in table 1, it can be seen that the 
optical and X-ray pulse profiles are approximately in phase 
with each other. To quantify this, the optical pulse profile 
was cross-correlated with the X-ray pulse profile. The result- 
ing peak in the cross-correlation function was fitted with a 
Gaussian to derive a phase shift of —0.06 ± 0.02 cycles (i.e. 
—0.39 ± 0.13s), where a negative phase shift implies that 
the X-ray pulse profile lags the optical pulse profile. This 
phase shift is only significant at the 3<r level, due to the low 
signal-to-noise and time resolution of the optical data, and 
additional data will be required in order to confirm that the 
phase shift is different from zero (discounting the unlikely 
situation in which the time delay is approximately equal to 
some multiple of the spin period). 

It should be noted that the morphology of the X-ray 
pulse profile in IE 1048.1-5937 does not appear to be en- 
ergy sensitive; the shapes of the 2-4 keV and 6-10 keV pulse 
profiles are virtually identical to the 2-10 keV pulse profile 
shown in figure [21 even though the 6-10 keV band is com- 
posed primarily of non-thermal photons whereas the 2-4 keV 
band is composed of both thermal and non-thermal photons 
(F. Gavriil, private communication). 

The modulation amplitude of the pulses presented in 
figure[2]can be measured using a peak-to-trough pulsed frac- 
tion, h p t, defined as follows: 
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frequency (Hz) 

Figure 3. Lomb-Scargle periodograms of IE 1048.1-5937 in the 
j'-band, obtained using the light curves from technique (ii) [sec- 
tion 12.21 . The dotted line shows the X-ray pulse frequency of 
0.1548479469 Hz given in table 1. 
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where F max and F m in are the maximum and minimum 
flux in the pulse profile, respectively. We find a value of 
hpt = 52 ± 15%. The peak-to-trough pulsed fraction defined 
in equation 1 effectively adds any noise present in the light 
curve to the true pulsed fraction, thereby tending to increase 
the resulting measurement. A more robust estimate is given 
by the root-mean-square (rms) pulsed fraction, h rrns , de- 
fined as follows: 



forms — _ 

y 



1 n 
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where n is the number of phase bins per cycle, yt is the 
number of counts in the i th phase bin, at is the error on 
yi and y is the mean number of counts in the cycle. As 
expected, measuring the optical pulsed fraction in this way 
gives a lower value of h rm s = 21 ± 7%. 

For comparison, the 2-10 keV X-ray rms pulsed frac- 
tion at the same epoch as the optical observations was 
hrms = 28.7 ± 0.5%. Since it isn't possible to measure the 
rms pulsed fraction directly from the RXTE pulse profile 
due to the uncertain background (see caption to figure [2]), 
we derived this value as follows. We first averaged the rms 
pulsed flux, defined as the product of the total flux and the 
rms pulsed fraction, measured with R XTE between 2007 
June 7 and June 12 (|Dib et alj|2008h . This value, which 
is background independent, was ide ntical to that me asured 
with Chandra on 2007 April 28 by iTam et ail (|2008T ). who 
also find a strong anticorrelation between total flux and rms 
pulsed fraction. Hence if the pulsed flux was the same for 
the RXTE and Chandra observations, we can be confident 
that the rms pulsed fraction was the same as well, and hence 
we have adopted the Chand ra rms pulsed frac tion from the 
2007 April 28 observation of ITam et all (|2008T ). 

3.2.2 Technique (ii) 

The second data reduction technique (section 12. 2[) can be 
used to provide a check on the reliability of the optical pulse 



profile derived using technique (i). Rather than adopting the 
X-ray ephemeris given in table 1, we can instead determine 
the pulse period directly from our optical data using a peri- 
odogram. Measuring the same period in the optical and the 
X-rays would prove without doubt that we have detected 
optical pulsations from IE 1048.1-5937. 

Figure B] shows the Lomb-Scargle periodograms 
l|Press fe Rvbickill 19891 ) for the 11 095 points in the i'-band 
light curves. The light curve was first corrected for trans- 
parency variations. The highest peak in the resulting pe- 
riodogram occurs at a frequency of 0.15484 ± 0.00004 Hz 
(6.458 ± 0.002 s), where the error is given by the width (a) 
of a Gaussian fit to the peak in the periodogram. This fre- 
quency is consistent to the fifth decimal place with the X- 
ray pulse frequency given in table 1, thereby confirming that 
we have indeed detected the X-ray pulsation of IE 1048.1- 
5937 in the optical. We further tested the robustness of our 
period detection by constructing 10 000 randomised light 
curves from the original light curves by randomly re-ordering 
the time-series. None of the resulting 10 000 periodograms 
showed a higher peak at 0.15484 Hz. To check for artifacts, 
we calculated the Lomb-Scargle periodograms of both the 
sky and the comparison star - neither showed any evidence 
for a periodicity at 0.15484 Hz, or at any other frequency. We 
also searched for evidence of optical bursts and/or longer- 
timescale periodicities in our 10 684 s light curve - nothing 
significant was found. 

Folding the 11 095 points of the technique (ii) light curve 
on the X-ray period given in table 1 results in the pulse 
profile shown by the dotted line in the top panel of figure [2] 
Note that the phasing of both the technique (i) and (ii) 
profiles in figure [2] can be directly compared, as all data 
have been folded using the zero point given in table 1. As 
one would expect, the technique (ii) light curve shown is 
in excellent agreement in terms of morphology, phasing and 
pulsed fraction (h rm s = 26 ± 8%) with the technique (i) 
light curve, lending additional confidence to our reduction 
and analysis techniques. 



4 DISCUSSION AND CONCLUSIONS 

The results presented in section [3] demonstrate conclusively 
that we have detected optical pulsations from IE 1048.1- 
5937 on the X-ray spin period. 

It is instructive to compare the optical light curve of 
IE 1048.1-5937 with that of the only other AXP to hav e 
been studied in this way - 4U 0142+61 (|Dhillon et al.ll2005t ). 
Both objects show optical pulsations on the X-ray spin pe- 
riod, which is 6.458 s in the case of IE 1048.1-5937 and 
8.688 s in 4U 0142+61. Both objects show optical pulsations 
with similar morphologies to their 2-10 keV X-ray light 
curves, with IE 1048.1-5937 exhibiting a single-humped pul- 
sation and 4U 0142+61 a double-humped pulsation. Both 
objects exhibit optical pulsations which are approximately 
in phase with the X-ray pulsations, with IE 1048.1-5937 
showing only marginal evidence for the optical leading the 
X-rays and 4U 0142+61 showing only marginal evidence 
for the optical lagging the X-rays. Even the optical pulsed 
fractions of the two objects are similar, with values of 
h pt = 52 ± 15% and h rms = 21 ± 7% in IE 1048.1-5937, 
and hpt = 58 ± 16% and h rm s = 29 ± 8% in 4U 0142+61. 
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The only major difference when comparing this study 
of IE 1048.1-5937 with Dhillon et al.'s (2005) study of 
4U 0142+61 is the ratio of the optical to X-ray pulsed frac- 
tion: in IE 1048.1-5937 it is approximately unity, whereas 
in 4U 0142+61 the optical pulsations had an rms pulsed 
fraction 5-7 times that of the X-rays. However, whereas 
the optical and X-ray pulsed fractions for IE 1048.1-5937 
were measured contemporaneously, those of 4U 0142+61 
were not. The X-ray pulsed fractions of 4U 0142+61 re- 
p orted bvlDhillon et all l|2005h were quoted from the work 
of iPatel et all (|2003l ). who obtained their Chandra data in 
2000, over two years prior to the opt i cal ob servations. We 
now know from the work of iDib et all (2007), however, that 
4U 0142+61 exhibited an increase in its pulsed light of 36% 
in the 2-10 keV band between 2002 and 2004. Hence it is 
possible that at least part of the discrepancy between the 
optical/X-ray pulsed fraction ratio in IE 1048.1-5937 and 
4U 0142+61 is due to variability in the latter so urce. This hy- 
pothe sis is further supported by the fact that lDhillon et all 
l|2005r ) meas ured a value of h v t = 58 ± 16% in September 
2002 whereas lKern fc Martini (|2002r i measured h pt = 27+8% 
in November 2001, although there are other pos sible rea- 
sons for this discrepancy (see iDhillon et all 120051 for a dis- 
cussion). It is also possible that the discrepancy between 
the optical/X-ray pulsed fraction ratio in IE 1048.1-5937 
and 4U 0142+61 is due to variability in IE 1048.1-5937: 
as discussed in section 13.2.11 the X-ray pulsed fraction of 
IE 1048.1-5937 is inversely proportional to the total 2-10 
keV flux. The optical pulsed fraction might not vary in the 
same way, implying that the optical/X-ray pulsed fraction 
ratio in IE 1048.1-5937 could be variable and we just hap- 
pened to observe it when it was unity. 

Viewed in isolation, the data on IE 1048.1-5937 pre- 
sented in this paper do not allow us to discriminate between 
the magnetar and fallback disc models, as the optical and X- 
ray pulsed fractions are approximately equal, although the 
tentative evidence we present for the optical pulses lead- 
ing the X-rays is irreconcilable with reprocessing. Arguably 
the most important result of this paper, however, is that 
it confirms the existence of pulsed optical light in a sec- 
ond AXP, IE 1048.1-5937, thereby demonstrating that the 
properties of the optical emission in 4U 0142+61 are not 
unique. The observed similarities between the optical and 
X-ray emission in both IE 1048.1-5937 and 4U 0142+61 in- 
dicate that closely related populations of particles, located 
in the same region of the magnetosphere, are probably re- 
sponsible for the emission. This paper has also highlighted 
the way forward for time-resolved optical studies of these 
incredibly faint objects, as the data presented in this paper 
would have been unobtainable had IE 1048.1-5937 been in 
a faint state. By targetting observations during bright X-ray 
states, it should be possible to study other AXPs, and pos- 
sibly also soft gamma repeaters, in the optical part of the 
spectrum, although this will still require access to sensitive, 
high-speed cameras like ULTRACAM on the world's largest 
telescopes. 
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